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Abstract 


Performance  of  a  grooved  heat  pipe  radiator  has  been  analyzed  to  determine  its 
operating  characteristics.  A  heat  transfer  analysis  for  the  heat  pipe  is  coupled  with 
the  analysis  of  the  fin  in  this  study.  The  effects  of  heat  pipe  operating  temperature 
on  heat  flux  and  fin  efficiency  are  investigated.  Finite  difference  expressions  are  used 
for  the  grid  system  of  the  heat  pipe  wall  and  the  fin. 

The  heat  transfer  coefficient  of  the  heat  pipe  radiator  was  determined  as  a 
function  of  the  operating  temperature  level.  Realistic  numerical  results  were  achieved 
and  it  was  shown  that  the  required  heat  flux  and  the  temperature  difference  between 
the  evaporator  and  condenser  both  increase  with  operating  temperature. 

Also,  the  heat  transfer  coefficient  increased  with  operating  temperature  up  to 
a  certain  point,  then  began  to  level-off  near  a  temperature  of  approximately  400  K. 
In  addition,  the  overall  efficiency  of  the  heat  pipe-fin  arrangement  decreased  as  the 
operating  temperature  level  increased. 
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THERMAL  ANALYSIS  OF  HEAT  PIPE  RADIATORS  WITH  A 
RECTANGULAR  GROOVE  WICK  STRUCTURE 


1.  INTRODUCTION 

1.1  Background 

Advanced  thermal  management  systems  are  necessary  for  future  military  appli¬ 
cations  in  space.  The  basic  components  of  these  thermal  systems  are  heat  aquisition, 
heat  transport,  and  heat  rejection  systems.  Current  power  levels  range  from  1  to  5 
kWg  and  heat  fluxes  are  less  than  1  W/cm^.  Future  power  levels  will  reach  50  -  100 
kWj  and  heat  fluxes  will  be  approximately  100  W/cm^. 

The  thermal  control  approach  employed  for  current  spacecraft  is  to  mount  the 
payload  on  the  spacecraft  structure.  Heat  is  transferred  by  heat  conduction  through 
the  structure  and  then  by  thermal  radiation  from  the  surface  of  the  structure  to 
space.  As  power  levels  and  thermal  control  requirements  increase  in  the  future, 
larger  heat  rejection  area  must  be  provided.  A  possible  solution  is  to  employ  a  heat 
pipe  radiator.  The  heat  pipe  radiator  may  be  designed  to  be  fixed  or  deployable 
with  two-sided  heat  rejection  capability.  The  deployable  radiator  may  also  have 
gimballing  mechanisms  to  lessen  radiation  heat  gains  from  the  sun  or  albedo. 

1.8  The  Heat  Pipe 

The  heat  pipe  was  developed  by  Grover  and  his  associates  for  use  in  spacecraft 
power  system[2,4],  and  the  first  analytical  work  was  carried  out  by  Cotter [2, 4].  The 
heat  pipe  is  a  closed  two-phase  thermal  device  which  operates  without  any  auxil¬ 
iary  external  pumping  power.  The  device  is  basically  composed  of  a  container,  a 
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capillay  structure,  and  a  working  fluid.  This  simple  thermal  device  is  able  to  trans¬ 
port  large  quantities  of  heat  with  a  very  small  temperature  drop  across  the  length. 
Furthermore,  it  is  versatile  in  shape  and  compact  in  size. 

The  heat  pipe  uses  phase  change  processes  to  accomplish  heat  transport.  A 
conventional  heat  pipe  is  depicted  in  Fig.1.1.  The  liquid  phase  of  the  working  fluid  is 
contained  in  the  capillary  structure  and  the  vapor  phase  placed  in  the  open  volume  of 
the  pipe.  As  heat  is  applied  to  one  surface  area,  the  so-called  evaporator,  the  liquid 
in  the  capillary  structure  starts  to  vaporize  and  thus  increases  the  vapor  pressure 
of  that  section  increases.  The  resulting  pressure  difference  between  the  ends  of  the 
pipe  causes  the  vapor  movement  carrying  the  latent  heat  of  vaporization  to  the 
other  end,  which  is  called  the  condenser.  As  heat  is  released  from  the  surface  of  the 
condenser,  vapor  is  condensed  in  the  capillary  structure.  This  condensate  returns  to 
the  evaporator  by  the  capillary  surface-tension  driving  force. 

By  choosing  a  suitable  working  fluid,  the  heat  pipe  can  operate  over  the  range 
of  temperatures:  cryogenic  (  0  -  150  K  ),  low  temperature  (  150-  750  K  ),  and  high 
temperature  (  750  -  3000  K  ).  The  operating  temperature  of  the  heat  pipe  must  be 
between  the  freezing  and  critical  temperatures  of  the  working  fluid,  because  success¬ 
ful  operation  needs  a  capillary  structure  saturated  with  liquid.  Heat  pipes  used  in 
this  study  for  the  radiator  use  grooves  for  the  capillary  structure.  The  working  fluid 
is  water  and  the  pipe  material  is  copper.  The  useful  operating  temperature  range  of 
water  heat  pipes  is  from  300  K  to  475  K.  For  this  temperature  range,  the  heat  trans¬ 
fer  through  the  liquid-saturated  capillary  structure  is  mainly  by  heat  conduction  and 
phase  change  takes  place  at  the  interface  between  the  liquid  and  vapor.  The  overall 
temperature  drop  thus  results  mainly  from  the  radial  temperature  drop  through  the 
pipe  wall  and  capillary  structure  in  the  evaporator  and  condenser  section. 


1-2 


l.S  Objectives  of  the  Study 

The  primary  purpose  of  this  study  has  been  to  analytically  investigate  normal 
operating  characteristics  of  a  water  heat  pipe  radiator  and  the  thermal  efficiency  of 
radiator  fins  as  a  function  of  the  radiation  parameter.  For  this  purpose,  a  mathe¬ 
matical  model  and  a  numerical  solution  technique  has  been  developed. 

Given  the  operating  heat  pipe  temperature  and  the  applied  heat  to  the  evap¬ 
orator,  temperatures  such  as  the  surface  temperatures  of  the  evaporator  and  con¬ 
denser  and  the  interface  temperatures  between  the  liquid  and  vapor  at  these  sections 
have  been  determined.  The  overall  heat  transfer  coefficient  has  also  been  found  as 
functions  of  the  operating  temperature  and  heat  input. 

Using  the  results  from  the  heat  pipe  analysis,  the  thermal  efficiency  of  the 
radiator  fins  has  been  investigated.  For  this  analysis,  it  is  assumed  that  the  radiator 
surface  is  a  gray  body,  which  has  the  same  value  for  the  emissivity  and  absorptivity. 

1.4  Literature  Survey 

Before  discussing  the  details  of  this  thesis,  a  review  of  the  literature  relating 
to  fin-radiator  studies  is  appropriate. 

Lieblien(1959)  studied  the  heat  transfer  characteristics  of  a  radiator  with  rec- 
tagular  fins.  It  was  cissumed  that  the  fin  had  a  uniform  temperature  at  the  edge, 
and  no  consideration  was  given  to  the  effects  of  fin-tube  mutual  irradiation. 

Sparrow,  Eckert,  and  Irvine(1961)  studied  the  general  case  of  a  radiator  with 
longitudinal  or  annular  fins,  which  included  mutual  irradiation  effects  between  the 
fin  and  tube,  which  was  reflected  in  their  governing  heat  transfer  relations. 

Annular  circular  fins  with  a  black  surface  attached  to  the  circular  tubes  was 
also  analyzed  by  Sparrow(1962)  with  consideration  of  mutual  irradiation  effects. 

Sparrow  and  Eckert(1962)  studied  the  effects  of  mutual  irradiation  occurring 
between  a  fin  and  a  base  surface  with  black  surfaces.  They  found  that  the  fin  heat 
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loss  Weis  significantly  reduced  by  the  presence  of  the  base  surfaces  in  the  range  of 
practical  operating  conditions.  Additionally,  the  beise  heat  surface  heat  loss  com¬ 
prised  an  important  part  of  the  total  heat  loss  from  the  given  system.  They  also 
formulated  the  equation  for  the  selective  gray  surfaces  and  arbitrary  irradiation  from 
external  sources,  and  analyzed  circular  black  tubes  with  longitudinal  fins  of  negligible 
thickness. 

Sarabia  and  Hitchcock(1966)  formulated  the  heat  transfer  equations  for  fin- 
tube  radiators  with  gray,  diffuse  emitting  and  absorbing  surfaces.  A  numerical 
method  was  also  developed  and  the  accuracy  of  the  correlation  between  the  gray 
and  black  surface  wa.s  established. 

The  total  weight  of  an  energy  conversion  system  is  very  important  for  space 
applications.  The  temperature  level  of  the  radiator  influences  the  total  weight  re¬ 
quired  to  dissipate  the  radiation  into  space.  For  design  studies  utilizing  finned-tube 
radiators,  the  results  of  references  14  and  3  can  be  used  to  determine  approximately 
the  optimum  operating  radiator  temperature  and  optimum  tube  radius  to  fin  length 
ratio. 

Rolf-Dieter  and  Walter(1962)  studied  the  minimum  configurations  of  tube- 
fin  radiators  with  gray  surfaces  and  the  governing  by  non-linear  integro- differential 
equation  which  can  be  solved  by  successive  approximations  using  Galerkin’s  method. 

From  the  literature  search,  it  is  evident  that  there  is  considerable  interest  in 
the  high  thermal  heat  transfer  radiator  for  space  applications.  Here,  the  study  of  the 
characteristics  of  the  heat  pipe-fin  radiator,  which  may  be  used  for  space  applications, 
has  been  carried  out  numerically. 
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II.  MATHEMATICAL  ANALYSIS 


Figure  2.1  depicts  a  grooved  schematic  of  the  heat  pipe  radiator  for  space 
applications.  Its  cross-sectional  view  is  shown  in  Fig. 2. 2,  which  presents  grooved 
heat  pipes  connected  by  rectangular  fins.  This  chapter  describes  a  mathematical 
analysis  under  the  following  assumptions: 

(a)  Heat  transfer  from  the  heat  pipe  condenser  surface  is  by  radiation  only. 

(b)  Heat  fluxes  into  the  evaporator  and  the  condenser  are  uniform. 

(c)  Heat  pipe  condenser  has  a  constant  temperature  distribution  along  the  heat 
pipe  axis. 

(d)  Heat  flow  in  the  heat  pipe  is  one  dimensional. 

(e)  Vapor  flow  in  the  heat  pipe  is  laminar  and  incompressible. 

(f)  Temperature  of  the  fin  base  is  constant  along  the  length  of  the  heat  pipe 
condenser. 

(g)  Heat  pipe  radiator  is  exposed  to  deep  space. 


2.1  Heat  Balance  Equation 


2.1.1  Heat  Pipe  Section  The  primary  heat  transfer  mechanisms  for  the  low- 
temperature  heat  pipe  are  (i)  heat  conduction  across  the  container  wall  and  the 
liquid-saturated  wick  at  the  evaporator  section,  (ii)  convective  transport  of  latent 
heat  by  vapor  from  the  evaporator  to  the  condenser,  and  (iii)  heat  conduction  across 
the  liquid-saturated  wick  and  the  container  wall  at  the  condenser  section.  The  heat 
conduction  in  the  pipe  wall  and  in  the  liquid-saturated  wick  can  be  described  by 
Fourier's  Law  ais  follows  : 


Q  = 


AT 

R 
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Heat  Pipe  Spacing(2L) 


Figure  2.1.  Multiple  Heat  Pipe  Radiator  in  the  Spacecraft 


Figure  2.2.  The  Cross-Sectional  View  of  the  Condenser  Section  of  Heat  Pipe 
Radiator 
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Evaporator  Section  Adiabatic  Section  Condenser  Section 


Figure  2.3.  Sketch  of  Heat  Flow  Path  through  a  Single  Heat  Pipe 

where  A  T  is  the  temperature  drop  through  the  wall  and  R  is  the  thermal  resistance, 
which  can  be  evaluated  from  expressions  for  the  cylinder  wall: 

M^2/ri) 

2w  ■  w  •  k 


where  k  is  the  thermal  conductivity  of  the  heat  pipe  wall  material  and  w  is  the 
length  of  heat  pipe  wall.  The  axial  heat  convection  by  vapor  can  be  described  by 
the  Clausius-Clapeyron  equation  which  relates  the  vapor  temperature  and  pressure 
as  follows  : 


Ti  -  Tz  = 


UPi  -  P2) 


Pv  '  hjg 

where  hjg  is  the  latent  heat  of  the  vaporization,  Ty  is  the  vapor  temperature  of 
working  fluid,  and  is  the  vapor  density  of  working  fluid. 
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Application  of  the  above  equations  to  each  section  in  a  series  of  heat  flow 
paths  in  a  conventional  heat  pipe  as  shown  in  Fig.  2.3,  with  the  temperatures  at  the 
different  locations  indicated,  yields 


•  Pipe  Wall  at  Evaporator  : 


T  —T  — 


2n  ■  We  ■  kp 


(2.1) 


where  kp  is  the  thermal  conductivity  of  the  heat  pipe  container,  and  Wj  is  the 
length  of  the  heat  pipe  evaporator. 

•  Liquid-Saturated  Wick  at  Evaporator  : 


T  —T  = 

pw^e  wv^e  — 


2n  ■  We  -  ke,e 


(2.2) 


where  ke,e  is  the  effective  thermal  conductivity  of  the  evaporator  capillary  struc¬ 
ture  (  See  Table  2.1). 


•  Axial  Vapor  Passage  : 


J  v,e  v.c  —  I  ^ 

Pv  "  "  V 


(2.3) 


•  Liquid-Saturated  Wick  at  Condenser  : 


T  —T 

WV.C  •*  pWtC 


JhAiIlALq 

2n  ■  Wc’  ke,c 


(2.4) 


where  ke,c  is  the  effective  thermal  conductivity  of  the  condenser  capillary  struc¬ 
ture  (  See  Table  2.1  ). 

•  Pipe  Wall  at  Condenser  : 


1  pw,c 


T  = 


2Tr  •  Wc  -  kp 


(2.5) 
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Table  2.1.  Expressions  of  Effective  Thermal  Conductivity  ke  for  Liquid-Saturated 
Wicks 


Rectangular  grooves 

kf  Expressions  ’ 

Evaporator 

ijw fWdkik-uj')  -j-  wki{Q.\hhw -j-  Wdki'j 

{wj  +  w){0.l85wjk^  +  Wdkt) 

Condenser 

_  wkt  Wfk^ 

W  -j-  Wf 

Where  h  =  liquid  thermal  conductivity 

k^  =  thermal  conductivity  of  wick  material 
w  =  groove  thickness 
Wf  —  groove  fin  thickness 
lud  =  groove  depth 


Combining  all  these  equations,  Eqs.(2.1)-(2.5),  results  in 
m  m  ,  /n(r./r„) 

p,e  p,c  ~  Vlo  I  ' 


27r  ■  We  ■  kp  2n  ■  We  ■  fce,e 

^  Tv[Pv,e  -  Pv,c)  ln{r,lr^)  /n(ro/r,) 
Pv  *  ^  ’  We  ' 

The  total  heat  transfer  of  the  heat  pipe  is  defined  as 


Q  =  AU„p[Tp,e  -  Tp,c) 


(2.6) 


(2,7) 


where  Tp.e  is  the  outer  surface  temperature  of  the  evaporator,  Tp,c  is  the  outer  surface 
temperature  of  the  condenser,  and  Upp  \s  the  overall  heat  transfer  coefficient  based 


on  an  arbitrary  area,  A.  Note  that  A  may  be  the  cross-sectional  area  of  the  pipe  Ap, 
the  surface  area  of  the  evaporator  A*,  the  surface  area  of  the  condenser  Ac,  or  any 
other  area.  By  nature  of  the  definition  for  Uhp,  it  follows  that 


AUhp  =  ApUpp^j,  =  AeUnP.e  =  AcUhp,c 


Choosing  the  cross-sectional  area  for  this  area,  Eq.(2.7)  becomes 


Q  =  ApUpp.piTp^e  -  Tp^c) 


where  Ap  is  Therefore,  from  Eqs.(2.7)  and  (2.8), 


(2.8) 


2  2 
0  °\  I  o 


rl  .  .r..  r:  .  ,  r,  T:rlT^{P^_^  - 


-Me)  + 


t,t  •  V 


Pv  ■  '  Q 


T'  T 


(2.9) 


2wc  •  ke^c  2u>c  •  kp  r, 

The  vapor  pressure  drop,  Pv^t  —  Pv,c  in  Eq.(2.9),  is  the  sum  of  the  average  vapor 
pressure  drops  in  the  evaporator,  the  adiabatic  section,  and  the  condenser. 


Pv,e  ~  Pv,c  —  ^Pv,e  +  ^Pv,a  +  ^Pv,c 


(2.10) 


These  components  of  the  vapor  pressure  drop  can  be  obtained  by  integrating  the 
following  equation  which  is  derived  from  the  axial  momentum  flux  law  [l]. 


dx 


=  -FvQ  -  Dy 


dQ^ 

dx 


(2.11) 


Here,  Fy  and  Dy  are  respectively  the  frictional  and  the  dynamic  pressure  coefficients 
for  the  vapor  flow,  which  are  dependent  upon  flow  conditions(  see  Table  2.2). 

When  the  vapor  flow  is  assumed  to  be  laminar  and  incompressible  with  neg¬ 
ligible  dynamic  effect  and  the  distribution  of  heat  flux  as  along  the  evaporator  and 
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Table  2.2.  Expressions  of  Vapor  Frictional  Coefficient  and  Dynamic  Coefficient 

D. 


For  circular  vapor 

core  cross  section 

Flow  conditions 

Fv 

Rcv  <  2300 

1.33 

M*,  <  0.2 

rlA^p^hfi, 

Alpyhjg 

Rcv  <  2300 
Mv>0.2 

)(1  + 

1.33 

Alpyh)^ 

Rcv  >  2300 
M^<0.2 

^  0.019/iv  -'^vQ  ^3/4 

AyjT^pXjhjg  .AyjhjgPy 

1 

Rcv  >  2300 

Mv>  0.2 

0.019mv  ^  2r„Q  3^  +  ~  ^ 

^  A,rlp^hjJ\A,h!gpJ  '  2 

1 

AlPuh}^ 

Rty  =  Reynolds  number  of  vapor  flow 
=  Mach  number 
Au  =  Vapor  flow  area(-ry) 
iL..  =  Vi.scositv  of  vapor 
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condenser  uniform,  the  average  pressure  drop  through  each  section  can  be  expressed 
as  follows: 

A  p 

D 

APv,a  =  F^WaQ 


Therefore, 


P.,.  -  P.,c  =  F,Q{^  +  w.  ^ 

6  6 


(2.12) 


Substituting  Eq,(2.r2)  into  Eq.(2.6)  gives 

'T  -  'T  ,  ln{r,fr^)  ^  T^F^{tL'e/&  +  ll'a  + 

p.c  —  J  p.e  ;  T  - ; - r 


‘2-  •  U'e  ■  kp  2TT  '  We  '  ^e.e 


Injrjrd  /^(rjr,) 

■  /-w  t  * 


Pv  •  h 


:s 


2t:  •  W'c  •  ke^c  27:  •  Wc  ■  kp 

2.1.2  Heat  Pipe  Condenser  Section  with  Fin 


(2.13) 


2. 1.2.1  Fin  Surface  The  basic  thermal  analysis  of  the  heat  pipe  con¬ 
denser  analysis  with  an  extended  fin  surface  is  the  conservation  of  energy  law.  The 
terminology  and  coordinate  system  for  the  longitudinal  fin  of  rectangular  profile  ra¬ 
diating  to  deep  space  are  shown  in  Fig.  2.4.  This  figure  also  shows  the  assumed 
direction  of  heat  flow  for  a  typical  fin  segment.  For  this  fin  segment,  symmetrical 
heat  transfer  from  both  (i.e., upper  and  lower)  surfaces  of  the  fin  is  assumed.  Under 
steady-state  conditions,  energy  conservation  for  the  differential  fin  element  dx  takes 
the  form  of  a  balance  between  net  conduction  and  net  radiation  effects. 


{dQcond)nrt  H"  {dQr  ad^net  —  0  (2. 14  ) 

The  fin  receives  heat  on  its  faces  from  both  heat  pipe  condenser  tubes.  Heat  enters 
into  the  fin  uniformly  at  its  base  where  x  =  0.  Referring  to  Fig. 2. 4,  the  heat  balance 
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Figure  2.4.  Heat  Flow  Path  in  a  Differential  Fin  Element 


equation  for  the  differential  element  dx  is 


Qc  +  /(x)  Wcdx  =  Oc  + 


dx  +  R{x)wcdx 
ax 


Fourier’s  law  for  one-dimensional  heat  conduction  gives 


dx 


-kj 


Wc  ■  t 


cPT 

dx^ 


(2.15) 


(2.16) 


where  kf  is  the  conductivity  of  the  fin,  which  is  made  from  the  same  material  ois  the 
heat  pipe,  and  lUc  is  the  length  of  heat  pipe  condenser.  Combining  Eqs.(2.15)  and 


(2.16)  yields 


I(x)wcdx  =  [—kj  •  Wc' t 


(FT 

dx'^ 


]dx  -f  R(x)wcdx 
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where 


=  A  ■ «/  = 


kf  ■  i 


The  boundary  conditions  are  for  the  governing  equation  as  follows: 


at  C  =  1, 


at  (  =  0, 


r*  =  1 


Integrating  Eq.(2.22)  once  gives 

—  no— rfcl 


(2.23) 


(2.24) 


(2.25) 


where 

A  =  0.6325/a 

and 

/■i  dT‘ 

'■  =  /o 

Finally,  the  temperature  at  the  new  axial  position  is  approximated  by  using  the 
Taylor  series  expansion, 

r(C  +  AC)  =  T-(0  +  ^|(AC  +  5^I<(AC)'  (2.26) 

where  the  derivatives  are  obtained  from  Eqs.(2.25)  and  (2.22). 

The  irradiation  of  the  fin  element  surface,  I(x),  consists  of  the  radiosity  received 
from  both  heat  pipe  condensers,  i.e.,  heat  pipe  condenser  1  and  condenser  2.  It  can 
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be  expressed  as 


I{x)wcdx  =  f  R[62)Fi^xXVcrod6\  +  f  R{62)F2^x'Wcf'od02  (2.27) 

JOl  J  $2 


In  dimensionless  form,  Eq.(2.27)  becomes 


/*(C)dC  =  /  R'{6,)F,^^rlde^  +  /  Rr{62)F2^irlde2  (2.28) 

JBi  J$2 


where 


!£ 

L 


2. 1.2. 2  Condenser  Base  Surface  Assuming  the  temperature  of  the  con¬ 
denser  surface  is  uniform,  an  expression  for  heat  loss  can  now  be  derived.  Considering 
the  condenser  surface  1  (  Fig. 2.2  ),  the  net  heat  loss  is  expressed  as  the  difference 
between  the  radiant  emission  and  the  incoming  energy  radiated  from  the  fin  and  the 
condenser  surface  2.  Therefore  the  heat  balance  equation  for  the  condenser  surface 
element  {wcrodO)  becomes 


q(0)wcrod0  -f  I{6)wcrod0  =  R{6)wcrod9 


(2.29) 


In  dimensionless  form,  Eq.(2.29)  becomes 


q\e)  =  R:{e)  -  r{e) 


(2.30) 


where 


q-{e)  = 


R^{e)  = 


m 
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and 
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In  dimensionless  form,  Eq.(2.36)  becomes 


ri0i)d0i  =  /  R-{x)Fi^^d9i  +  I  R\e2)F2^idd2  (2.37) 

J  ^  J$2 


where 


R'ix)  = 


2.2  Heat  Transfer 

2.2.1  Fin  Surface  The  fin  heat  transfer  rate  can  be  expressed  by 


^  .  dT, 

Q f  —  ^  |i=0 


(2.38) 


The  fin  efficiency,  r//,,,,  is  defined  the  ratio  of  the  heat  transferred  by  the  fin  to 
the  heat  that  could  be  dissipated  by  an  isothermal  black  fin  of  the  same  dimensions, 
and  can  be  expressed  as  follows: 


from  Eq.(2.38) 


Vfin  — 


Qf 


Vfin  — 


kf  •  Wc  ■  t  dT 
crT*^  -Wc’Ldx 


Vfi 


1  dT* 
'01  dC 


lc=o 


(2.39) 


where 


anrf  = 

2.2.2  Total  Surface  The  total  heat  transfer  rate  can  be  expressed  by 


Q total  =  q{0)worod0  +  Qj] 


(2.40) 
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The  overall  fin-condenser  efficiency,  •qejj  ,  is  defined  the  ratio  of  the  heat  trans¬ 
ferred  from  the  heat  pipe  condenser  and  the  fin  to  the  heat  that  could  be  dissipated 
by  a  fiat,  black  surface  of  the  same  total  length  the  heat  pipe  condenser  and  the 
fin.  Therefore, 

..  .  (2.41) 


i[cT^AL  +  Jrl  - 


Combining  Eqs.(2.40)  and  (2.41)  gives 


n^JJ  = 


fe{q{0)wcro)d9  -f  Qf 
crT*^,Wc[L  +  Jrl  -  t^] 


where 


5  =  1  +  ,/r*2  - 


Substituting  Eq.(2.33)  into  Eq.(2.42)  yields 


(2.42) 


where 


Vtff  =  Ctl-  -  Omin  -  Dt]  + 


1  +  Jr-J  - 


9min  =  tan  !  I  . J 


(2.43) 


Dt  =  r  r{9)d9 
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III.  COMPUTATIONAL  PROCEDURE 


3.1  Summary  of  Equations  and  Finite  Difference  Forms 
Eq.(2.12)  : 

b  b 


Eq.(2.3)  ; 


rp  T  —  ~  n 

Jf,e  ~  J-v.c  —  7  77  V 

Pv  ■  i^Jg  ■  Q 


Eq.(2.2)  ; 


Jq.(2.9) 


T  =T 

■^pw^e  —  r  ^ 


27r  * 


Eq.(2.1) 


^  ^  tn(rolri)  ^ 

Tp,t  =  Ppuj.e  +  o_  . . — ~Q 


27r  ■  We  ■  kp 


Eq.(2.4) 


rr  rr.  ln{rJrA  ^ 

'^pWyC  —  '^WV.C  ^  W 


2it  ■  Wc-  ke,c 


Eq.(2.5) 


^  ^  ln{rjr,)  ^ 

Ap,C  —  J-pW.C  o  7  V 

zn  •  xvc '  kp 


Eq.(2.6) 


217  ■  We  •  kp  2i7  ■  We  ■  ke,e 
Tv{Pv,e  -  Pv.c)  /n(r,/r„)  /n(ro/r,) 

»  4^  f  ' 


Pv-  hfg-Q  2ir  -Wc-  ke,c  2ir  ■  Wc  ■  kp 


Uhp,p  = 


2w.  ■  k, 


■ln{rolri)  + 


2We  •  ke,e 


ln{ri/r^)  + 


i7rlT,iP,,e  -  P,c) 

Pv  ■  hfg  ■  Q 
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From  Eq.(2.19) 


Eq.(2.27)  : 


B 

kf  •  t 


r*(C)c/c=  /  R'{6,)F^^,Tide,+  I  R’ie2)F2. 

Jei  JSi 


>xro(i^2 


In  a  finite  difference  form  ,  this  becomes 


where 


*/2  jr/2 

=  E  R\0,)Fr^,Ar^  +  E  R'{02)F2^xAr2 


^  ,A01  .  .A^2 

^ri=r,  — ,  ^^2  =  r,— 


0min  —  tan  (■ 


Eq.(2.21) 


-  /?/-(C) 


Eq.(2.24)  :  reduces  into  a  finite  difference  equation  as 

jT-i*  C  riT* 

—  =  -A{{T-'  -  Te^)  +  5{/;  -  EriO^AOl'" 
“C  0 


where 


A  =  0.6325\/^ 

/r  =  Enc)[^]Ac 
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Eq.(2.25)  : 

r(c  +  AC)  =  r(c)  +  ^|(Ac  +  i^lc(AC)" 

Eq.(2.20)  : 

R’{0  =  ejr^+[l-ej]r{0 


Eq.(2.36) 

ri0i)dei  =  f  R^iOFr^cdOi  +  /  R^{02)F2^^de2 

J  $1  J  $2 

Converting  into  the  finite  difference  equation, 


r(9.)  =  E^'(0^’.-< 
0 


ir/2 

+  ^  R\e2)F2^, 


e 


min 


A02 


Eq.(2.31)  ; 

fl-(«)  =  <p  +  (i -(,)/■(«) 

From  Eq.(2.42) 

IT 

’!■//  =  c.)-  -  -  E  r(«)A«)  + 

^min 

where 

r  - 
‘  B 

and 

5  =  1  +  v^r*2  - 

S.2  Initial  Values 

(3.3) 


(3.4) 


Because  of  symmetry,  the  midpoint  of  the  fin  is  considered  to  be  adiabatic. 
Therefore,  the  energy  balance  for  this  differential  control  volume  is 


(.jcT^Wcdx  =  2  R{0)F\-,EWc‘rod0 


(3.5) 
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or,  in  dimesionless  form, 


=  2r:  r  R*{e)F,^Ede  (3.6) 

The  initial  condenser  surface  radiation  is  eissumed  to  be  equal  to  the  condenser 
surface  emission  : 

R{0)wcrodd  =  tpaT^^^WcXodO  (3-7) 

In  dimensionless  form,  Eq.(2.21)  becomes 

R'{B)  =  tp  (3.8) 

Substituting  Eq.(2.24)  into  Eq.(3.1)  yields 

=  2r;cp  Fr^EdO  (3.9) 

^min 

Expressing  Eq.(3.9)  in  finite  difference  form  yields 

«m.n 

or 

2 

r|=1.189(e,r;— (3.11) 

Even  if  the  irradiation  on  the  fin  midpoint,  as  expressed  in  Eq.(3.3),  is  small, the 
lowest  value  that  the  fin  midpoint  temperature  can  have  must  still  be  larger  than 
that  for  the  case  of  fins  with  no  irradiation.  The  solution  to  the  simple  fin  problem  is 
included  in  Appendix  D.  A  fifth  order  polynomial  expression  for  the  simple  fin  edge 
temperature(  fin  midpoint  )  as  a  function  of  the  radiation  parameter,  (3,  is  given  in 
Eq.(D.3),  Appendix  D. 
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The  initial  midpoint  temperature  assumed  in  this  solution  was  the  one  cal¬ 
culated  by  Eq.(3.3)  or  Eq.(D.3),  whichever  was  greater.  The  initial  temperature 
gradient  along  the  fin  was  computed  with  the  following  equations  : 


and 


dT* 

dC 


\i  =  -A[T’^  - 


(3.12) 


r(c  +  AC)l.  =  r-(c)  +  ^|.<iC 


(3.13) 


Equation(3.2)  shows  that  the  temperature  gradient  depends  on  the  value  of 
It  which  in  turn  depends  on  the  product  of  the  irradiation  and  the  temperature 
gradient.  The  initial  values  of  the  temperature  gradient,  Eq.(3.12),  are  too  large 
near  the  fin  base.  Therefore,  in  order  to  attenuate  this  initial  error,  the  following 
linear  reduction  factor,  Rf,  was  used  for  the  first  10  iterations  and  applied  to  the  fin 
irradiation  : 


r(C)  =  /*(c)i?/  (3.14) 

where 


Rf  =  — 

^  10 

for 

n  <  10 

o 

II 

for 

n  >  10 

In  order  to  preclude  obtaining  results  that  include  the  reduction  factor,  the  crite¬ 
ria  for  ending  a  solution  included  the  requirement  that  at  least  15  iterations  be 
performed  in  all  cases. 


3.3  Method  of  Calculation 

Step  1  :  Choose  an  appropriate  operating  temperature  level. 


Step  2  :  Assume  T’p,c  is  approximately  the  same  as  the  given  operating 

temperature. 

Step  3  :  As  a  first  guess,  estimate  the  required  heat  input  using  Eq.  (2.41). 

Step  4  ■  Solve  Eq.  (2.12)  to  determine  the  vapor  pressure  drop. 

Step  5  :  Solve  Eq.  (2.3)  to  determine  the  temperature  difference  between 

evaporator  and  condenser  surfaces. 

Step  6  :  Solve  Eqs.  (2.1),  (2.2),  (2.4),  (2.5),  and  (2.6)  for  temperatures 

along  the  heat  pipe. 

Step  7  ;  Solve  Eq.  (2.9)  to  determine  the  overall  heat  pipe  heat  transfer 
coefficient. 

Step  8  :  Calculate  the  view  factors  and  the  initial  values. 

Step  9  :  Solve  Eq.  (3.1)  with  ^  ranging  from  0.0  to  1.0. 

Step  10  :  Solve  Eqs.  (2.22),  (2.25),  and  (2.26)  sequentially  at  each  value  of 
C  from  0.0  to  1.0.  The  fin  midpoint  temperature  used  in  these  equations  is  obtained 
from  the  initial  values  or  from  the  previous  iteration.  There  is,  therefore,  an  ambigu¬ 
ity  between  the  computed  and  the  assumed  midpoint  temperatures.  To  remove  this 
ambiguity,  the  solution  to  this  set  of  equations  was  repeated  using  a  fin  midpoint 
temperature  equal  to  the  average  between  the  computed  and  the  asssumed  values. 

Step  11  :  Solve  Eqs.  (2.21),  (2.32),  and  (3.3)  with  ranging  from  0.0  to  1.0 
and  6  from  9min  to  7r/2.0 

Step  12  :  Test  the  solution  for  the  requirement  that  the  temperature  gradient 
at  the  fin  midpoint  be  1.5  percent  of  the  fin  root  temperature  gradient  or  less.  If  the 
criteria  was  not  met,  the  iteration  is  continued  by  returning  to  Step  6.  The  use  of  a 
similar  criteria  in  the  solution  of  the  simple  fin  radiator  resulted  in  an  error  in  the 
fin  efficiency  of  about  one  percent. 
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Step  13  :  With  updated  value  for  T’p,c,  caculate  the  new  estimate  for  the 

required  heat  input. 

Step  14  :  Calculate  the  relative  error  using  the  new  and  previously  obtained 
heat  input  estimates.  If  this  error  is  less  than  a  certain  tolerance,  the  solution  is 
converged  and  the  calculated  heat  input  is  considered  the  correct  match  for  the  given 
operating  temperature  level.  If  the  solution  does  not  converge,  update  the  value  for 
the  heat  input  and  repeat  the  iteration  sequence  starting  with  Step  4- 
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IV.  RESULTS 


The  numerical  results  in  this  study  were  obtained  by  considering  different  heat 
pipe  operating  temperatures  and  then  calculating  the  corresponding  heat  rejection 
requirements.  The  finite-difference  scheme  used  a  system  of  100  increments  along  the 
fin  and  30  increments  around  one-fourth  of  the  heatpipe  condenser  circumference, 
cissuming  symmetry. 

4-i  Heat  Transfer  Coefficient  of  Heat  Pipe  Radiator 

Figure  4.1  shows  the  variation  of  the  radiation  parameter  with  respect  to  the 
temperature  of  the  heat  pipe  condenser  surface.  Fig.  4.2  shows  the  total  heat  transfer 
rate  of  the  heat  pipe  as  a  function  of  the  operating  temperature.  Fig.  4.3  shows 
the  temperature  difference  between  the  evaporator  and  the  condenser  surfaces  of 
the  heat  pipe  for  various  operating  temperature  levels.  By  using  Eq.(2.9)  the  heat 
transfer  coefficient  of  the  heat  pipe  radiator  was  also  determined  as  a  function  of  the 
operating  temperature  level.  The  numerical  results  are  tabulated  in  Tables  C.l,  C.2 
and  C.3  in  Appendix  C. 

4-1 -1  Operating  Temperature  Effect:  As  shown  in  Fig.  4.2,  as  the  operating 
temperature  increases  so  does  the  required  heat  transfer  level.  The  temperature 
difference  also  follows  a  similar  trend  (  see  Fig.  4.3  ),  but  the  rate  of  increase  was 
slightly  slower.  The  temperature  effect  on  the  heat  transfer  coefficient  is  given  in 
Fig.  4.4.  It  is  interesting  to  note  that  the  increase  in  heat  transfer  coefficient  begins 
to  level-off  at  approximately  400  K,  reaching  a  value  near  62,000.  Also,  as  the 
operating  temperature  is  raised,  the  radiation  parameter  increases,  thus  causing  the 
fin  efficiency  to  derease,  see  Fig.  4.1. 
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4- 1-2  Heat  Input  Effect:  As  the  heat  input  increases  with  operating  temper¬ 
ature,  the  overall  temperature  difference  between  the  evaporator  and  the  condenser 
also  increases,  see  Fig.  4.3. 

4-2  Black  Surfaces  with  Negligible  Fin  thickness 

Figure  4.5  shows  the  overall  heatpipe-fin  efficiency,  Tjeffi  ^  function  of  the 
radiation  parameter  (/5i).  The  numerical  results  are  included  in  Table  C.4  of  Ap¬ 
pendix  C,  and  compared  well  with  the  results  of  reference  14  in  Figs.  4.9,  4.10,  and 
4.11  [14]. 

4-3  Heat  Pipe-Fin  Radiator  Efficiency 

As  shown  in  Figs.  4.5,  4.6,  4.7,  and  4.8,  the  heat  pipe-fin  radiator  efficiency  de¬ 
creases  as  the  radiation  parameter  increases.  The  radiation  parameter  is  a  function 
of  condenser  temperature  which  is  determined  from  the  given  the  operating  temper¬ 
ature.  For  the  given  heat  input,  the  radiation  parameter  increaises  with  increasing 
operating  temperature,  and  this  has  two  opposing  effects  on  the  heat  pipe-fin  effi¬ 
ciency:  (i)  the  heat  is  dissipated  more  in  the  section  near  the  fin  base  than  in  the 
remaining  section  of  the  fin.  (ii)  Radiosity  from  the  heat  pipe  condenser  also  in¬ 
creases,  which  means  the  irradiation  on  the  fin  increases,  resulting  in  a  higher  fin 
temperature. 
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Figure  4.1.  Relationship  between  Radiation  Parameter  and  Operating 
Temperature 
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Figure  4.2.  Total  heat  transfer  rate  of  heat  pipe  vs.  Operating  Temperature 
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Figure  4.4.  Overall  heat  transfer  coefficient  of  heat  pipe  vs.  Operating  Temperature 
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Figure  4.7.  Simple  Fin  Efficiency  from  Reference  14 


Figure  4.8.  Condenser  Surfao 


4- 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  mathematical  analysis  for  a  heat  pipe  radiator  has  been  carried  out  to  in¬ 
vestigate  fin  efficiency  and  heat  rejection  requirements  as  a  function  of  the  heat  pipe 
operating  temperature.  Grooved  copper  heat  pipes  with  water  ais  a  working  fluid, 
and  copper  fin,  were  considered  for  a  given  heat  pipe  operating  temperature,  the 
surface  temperature  of  the  evaporator  and  the  condenser  sections  were  calculated, 
along  with  the  required  heat  rejection  level. 

From  the  evaluated  condenser  surface  temperature,  the  radiation  parameter 
and  in  turn,  the  fin  efficiency  were  found.  It  was  found  that  the  overall  temperature 
difference  between  the  evaporator  and  the  condenser  increases  as  the  operating  tem¬ 
perature  increases.  As  the  condenser  temperature  increases,  the  radiation  parameter 
increases  while  the  fin  efficiency  decreases. 

It  is  recommended  to  perform  a  more  involved  parametric  study  by  using  dif¬ 
ferent  dimensions  of  heat  pipes  and  fins.  Various  vapor  flow  conditions,  such  as 
incompressible  turbulent,  compressible  laminar,  and  compressible  turbulent,  may  be 
included.  Furthermore,  different  wick  structures  for  the  heat  pipe  may  be  incor¬ 
porated,  including  screens,  metal  fibers,and  metal  powders.  The  existing  computer 
program  was  written  in  a  general  fashion  to  facilitate  condition  changes  for  an  effec¬ 
tive  parametric  study. 


5-1 


Appendix  A.  FORTRAN  PROGRAM 
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* 

*  FORTRAN  PROGRAM  LISTING 

* 


*  VARIABLE  DEFINITIONS 


* 

* 

♦ 

* 

* 

♦ 

* 


* 

* 


♦ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


CKL  THERMAL  CONDUCTIVITY  OF  LIQUID  AT  CONDENSER  [  W  /  M*K  ] 

CKP  THERMAL  CONDUCTIVITY  OF  HEATPIPE  MATERIAL 

AT  CONDENSER  [  W  /  M*K  ] 

CKW  THERMAL  CONDUCTIVITY  OF  WICK  MATERIAL  AT  CONDENSER 

[  H  /  M*K  ] 

EKL  THERMAL  CONDUCTIVITY  OF  LIQUID  AT  EVAPORATOR  [  W  /  M*K  ] 

EKP  THERMAL  CONDUCTIVITY  OF  HEATPIPE  MATERIAL 

AT  EVAPORATOR  [  W  /  M*K  ] 

EKW  THERMAL  CONDUCTIVITY  OF  HICK  MATERIAL  AT  EVAPORATOR 

C  W  /  M*K  ] 

ETG  INITIAL  GUESS  VALUE  OF  OVERALL  EFFICIENCY 

FLEN  FIN  LENGTH  [  M  ] 

PKPC  THERMAL  CONDUCTIVITY  OF  HEAT  PIPE  CONDENSER  [  W  /  M*K  ] 
RI  HEATPIPE  INSIDE  RADIUS  [  M  ] 

RRO  HEATPIPE  OUTSIDE  RADIUS  [  M  ] 

RV  HEATPIPE  VAPOR  CORE  RADIUS  [  M  ] 

W  GROOVE  WIDTH  [  M  ] 

WA  LENGTH  OF  HEAT  PIPE  ADIABATIC  SECTION  [  M  ] 

WC  LENGTH  OF  HEAT  PIPE  CONDENSER  [  M  ] 

WD  GROOVE  DEPTH  [  M  ] 

WE  LENGTH  OF  HEAT  PIPE  EVAPORATOR  [  M  ] 

WF  LENGTH  OF  GROOVE  FIN 

NG  NUMBER  OF  RECTANGULAR  GROOVES 

TS  TEMPERATURE  OF  SPACE  ENVIRONMENT  [  K  ] 

Q  TOTAL  REQUIRED  HEAT  REJECTION  [  W  ] 

SIGMA  STEFAN-BOLTZMANN  CONSTANT  ;  5.670E-08  [  W  /  M**2  K**4  ] 
NHP  NUMBER  OF  HEAT  PIPE  IN  MULTIPLE  HEAT  PIPE  RADIATORS 

QREQ  HEAT  TRANSFER  RATE  IN  SINGLE  HEAT  PIPE  [  W  ] 

TPC  CONDENSER  PIPE  WALL  TEMPERATURE  [  K  ] 

PKEE  EFFECTIVE  THERMAL  CONDUCTIVITY  OF  LIQUID-SATURATED  WICK 

AT  EVAPORATOR  [  W  /  M*K  ] 

TPWE  EVAPORATOR  TEMPERATURE  AT  PIPE-HICK  INTERFACE  [  K  ] 

TWVE  EVAOPRATOR  TEMPERATURE  AT  HICK-VAPOR  INTERFACE  [  K  ] 

TVE  VAPOR  TEMPERATURE  IN  EVAPORATOR  [  K  ] 

VISCOSV  VAPOR  DYNAMIC  VISCOSITY  [  KG  /  S*M  ] 

VDENS  VAPOR  DENSITY  [  KG  /  M**3  ] 

TV  VAPOR  TEMPERATURE  C  K  ] 
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* 

* 

* 

* 

* 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

♦ 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 

* 


HFG  LATENT  HEAT  OF  VAPORIZATION  [  J  /  KG  ] 

AREAV  VAPOR  CORE  CROSS-SECTIONAL  AREA  [  M**2  ] 

AREAP  HEAT  PIPE  CROSS-SECTIONAL  AREA  [  M**2  ] 

AREAC  CONDENSER  SURFACE  AREA  OF  HEAT  PIPE 
FV  FRICTIONAL  COEFFICEINT  FOR  VAPOR  FLOW 

PDELTA  VAPOR  PRESSURE  DROP  [  N  /  M**2  ] 

TVC  VAPOR  TEMPERATURE  IN  CONDENSER  [  K  ] 

TWVC  CONDENSER  TEMPERATURE  AT  HICK-VAPOR  INTERFACE  [  K  ] 

PKEC  EFFECTIVE  THERMAL  CONDUCTIVITY  OF  LIQUID-SATURATED  WICK 

AT  CONDENSER  [  W  /  M*K  ] 

TPHC  CONDENSER  TEMPERATURE  AT  HICK-PIPE  INTERFACE  [  K  ] 

TPCNEW  NEW  CONDENSER  PIPE  HALL  TEMPERATURE  [  K  ] 

EPSIP  EKISSIVITY  OF  HEAT  PIPE  WALL  SURFACE  [  DIMENSIONLESS  ] 
THICK  FIN  THICKNESS  [  M  ] 

BETAl  RADIATION  PARAMETER  [  DIMENSIONLESS  ] 

UHPP  HEAT  TRANSFER  COEFFICIENT  BASED  ON  ’AREAP’  [  W/M**2-K  ] 

EPSIF  EMISSIVITY  OF  FIN  SURFACE  [  DIMENSIONLESS  ] 

NFIN  NUMBER  OF  FIN  NODES 

NPIPE  NUMBER  OF  HEAT  PIPE  NODES 

TSTAR  RATIO  OF  FIN  NODE  TEMPERATURE  TO  CONDENSER  HEAT  PIPE 
WALL  TEMPERATURE 

RSTAR  TATIO  OF  HEAT  PIPE  OUTSIDE  RADIUS  TO  HALF  LENGTH  OF  FIN 
RHOP  REFLECTIVITY  OF  HEAT  PIPE  WALL  SURFACE 
RHOF  REFLECTIVITY  OF  FIN  SURFACE 

VFPF  VIEW  FACTOR  FROM  FIN  SEGMENT  TO  HEAT  PIPE  SEGMENT 

VFPP  VIEW  FACTOR  FROM  HEAT  PIPE  SEGMENT  TO  HEAT  PIPE  SEGMENT 

TF  FIN  NODE  TEMPERATURE 

FIRRAD  IRRADIATION  ON  THE  FIN  AT  FIN  NODES 

FRAD  RADIOSITY  FROM  THE  FIN 

PRAD  RADIOSITY  FROM  THE  HEAT  PIPE 

TEMISMP  TOTAL  EMISSION  FROM  THE  HEAT  PIPES  TO  THE  FIN  MIDPOINT 
TEDGE  MIDPOINT  TEMPERATURE  OF  FIN 

RATIO  RATIO  OF  TEMPERRATURE  GRADIENT  AT  THE  MID  FIN  TO  THE 
ROOT  TEMPERATURE  GRADIENT 
RF  REDUCTION  FACTOR 

T2IRRAD  TOTAL  IRRADIATION  ON  A  SEGMENT  OF  ONE  HEAT  PIPE  FROM  THE 
OTHER  HEAT  PIPE 

TIRRAD  TOTAL  IRRADIATION  ON  A  HEAT  PIPE  SEGMENT  FROM  THE  FIN 
QTIRRAD  TOTAL  IRRADIATION  ON  THE  HEAT  PIPE  PER  UNIT  RADIUS 
ETAEFF  EFFICIENCY  OF  HEATPIPE-FIN  RADIATOR 

ETAEFF  EFFICIENCY  OF  FIN  ONLY 
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*  DIMENSIONS  AND  FORMATS 

* 

REAL  RRO , TPC , EKL , EKP , TV , VDENS , VISCOSV , HFG . TWVC , CKL . CKP , PKPC , WF 
REAL  RI ,  RV ,  WE ,  WA ,  WC ,  W ,  WD ,  TS ,  TPE ,  Q ,  FLEN ,  THI CK ,  EPSIP ,  EPS  I F 
REAL  BETA,BETA1,UHPP,PDELTA 

REAL  VSINEOl)  ,VC0SIN(31)  ,VFPF(101 ,31)  ,VFPP(31,31)  ,TF(101) 

REAL  GRADT(lOl) ,GRADT2(101) ,FIRRAD(101) ,FRAD(101) ,PRAD(101) 
INTEGER  PIPENO,FINNO,NG,NHP.NFIN,NPIPE 

* 

OPEN (UNIT=3 , FILE=  *  PIP 1 . DAT ’ , STATUS= » OLD ' ) 

OPEN (UNIT=6 , FILE= ’ PIP2 . DAT ’ , STATUS=  *  NEW » ) 

READ(3.91)  RRO.RV 
READ(3.90)  WE.WA.WC 
READ(3,91)  W,WD 
READ (3. 92)  NG.NHP 
READ(3,90)  TS.TV.ETG 
READ (3, 91)  FLEN, THICK 
READ(3,91)  EPSIP, EPSIF 
READ(3,92)  NFIN.NPIPE 

90  F0RMAT(3F12.6) 

91  F0RMAT(2F12.6) 

92  F0RMAT(2I5) 

*  CONSTANTS 

DC  4^  He  %  41 *****  >t> «  4<  *  >•>«**  4:  >X  *««****«* 


TPCG  =  TV 
PI  =  3.141593 
HFLEN  =  FLEN/2. 

HTHICK  =  THICK/2. 

SIGMA  =  5.670E-8 
RI  =  RV  +  WD 

WF  =  (  PI*(RI+RV)/FLOAT(NG))  -  W 

QTTl  =  (SIGMA*TPCG**4.*(HFLEN+SQRT(RR0**2.-HTHICK**2.))*WC) 
QREQ  =  4.*ETG4qTTl 
PRINT*,  QREQ 
799  CONSTl  =  RRO/RI 
C0NST2  =  RI/RV 

CONSTA  =  (  LOG(CONSTl)  *  QREQ  )  /  (  2.0*PI  ) 

CONSTB  =  (  L0G(C0NST2)  *  QREQ  )  /  (  2.0+PI  ) 

AREAV  =  PI*RV**2 
AREAP  =  FI*RR0**2 
AREAC  =  2.0*PI*RR0*WC 
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* 

*  PART  1  :  HEAT  PIPE  CONDITIONS 

* 

* 

*  . CALCULATE  PROPERTIES  t  PRESSURE  DROP 

*  AT  VAPOR  TEMPERATURE 

CALL  STEAMDENCTV.VDENS) 

CALL  STEAMUG(TV,VISCOSV) 

CALL  STEAMHFGCTV.HFG) 

*  . 

FV  =  (  8.0+VISCOSV  )  /  (  HFG*VDENS*AREAV*RV**2  ) 
PDELTA  =  FV*QREQ*(  HE/6.0  +  HA  +  HC/6.0  ) 

♦  . CALCULATE  TEMPERATURE  AT  EVAPORATOR 

TVC  =  TV 

TVE  =  TVC  +  (  TV*PDELTA  )  /  (  VDENS*HFG  ) 

THVE  =  TVE 

♦  . CALCULATE  CONDUCTIVITY  t  TEMPERATURE  AT  EVAPORATOR 

CALL  HATERCTHVE.EKL) 

CALL  COPPER (THVE, EKP) 

EKH  =  EKP 

PKH  =  0.185*HF*EKH  +  HD*EKL 

PKEE  =  (HD*HF*EKL*EKH  +  H*EKL*PKH)  /  ((H+HF)*PKH) 
TPHE  =  THVE  +  CONSTB  /  (  PKEE*HE  ) 

CALL  COPPER (TPHE, EKP) 

EKH  =  EKP 

TPE  =  TPHE  +  CONSTA  /  (  EKP+HE  ) 

♦  . CALCULATE  TEMPERATURE  AT  CONDENSER 

TVC  =  TV 

THVC  =  TVC 

CALL  HATER(THVC,CKL) 

CALL  COPPER(THVC,CKP) 

CKH  =  CKP 

PKEC  *  (  H*CKL  +  HF*CKH  )/(  H  +  HF  ) 

TPHC  =  THVC  -  CONSTB  /  (  PKEC*HC  ) 

TPC  =  TPHC  -  CONSTA/  (  CKP*HC  ) 

CALL  COPPER(TPC,CKP) 
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CALL  COPPERCTPC.PKPC) 


UHPP  =  (  QREQ  )  /  ((  TPE  -  TPC  )*  AREAP  ) 

BETAl  =  (  SIGMA* (TPC**3)*(HFLEN**2)*EPSIP  )/(  PKPC*HTHICK  ) 
BETA  =  BETA1*EPSIF 


WRITE(6,20)  RRO.RI.RV 

20  FORMAT ( IX, ’RRO  =’.F8.2,  SX.’Ri  =’ ,F8.2,5X, 'Rv  =',F8.2) 

WRITE(6,21)  FLEN, THICK 

21  FORMAT (IX, ’FLEN  =’,F8.2,  5X, ’THICK  =’,F8.2) 

WRITE(6,22)  WE,WA,WC 

22  FORMATClX, ’WE  =’,F8.2,  5X,’WA  =’ ,F8.2,5X, ’WC  =’ ,F8.2) 

WRITE(6,23)  W,HD,NG 

23  FORMATClX, ’W  =’,F8.2,  SX.’WD  *’ ,F8.2,5X, ’NG  =',I3) 

WRITE (6, 24)  EPSIF,EPSIP 

24  FORMAT ( IX, ’EPSIF  *’,F8.2,  SX.’EPSIP  =’,F8.2) 

WRITE(6,25)  NFIN,NPIPE.NHP 

25  FORMATClX, ’NFIN  =’,I3,  SX.’NPIPE  =’ ,I3,5X, ’NHP  =’,I3) 

WRITE(6,26)  TPE,TS,Q 

26  FORMATClX, ’TPE  =’,F8.2,  5X,’TS  =’ .F8.2.5X, ’Q  =’,F8.2) 

WRITE(6,38) 

WRITE ( 6 , 39 )  TPE , TPWE , TVE , TV , TPWC , TVC , TPC 
WRITE (6, 40)  TPE, TPC 
WRITE(6,41)  BETAl, UHPP, PDELTA 
WRITE (6, 42) 

WRITE(6,43)  TPE, EKP,EKL, TPWE, TVE 
WRITEC6 ,45)  TEMPDIFF , BETAl .UHPP ,PKPC 
WRITE(6,44) 

WRITE(6,46)  TV, HFG, TVC, TPWC, TPC 

38  F0RMAT(8X, ’TPE’ ,6X, ’TPWE’ ,7X, ’TVE’ ,8X, ’TV’ ,6X, ’TPWC’ ,7X, ’TVC’ , 

1  7X,’TPC’5X) 

39  FORMATClX,  7(2X,F8.2)) 

40  FORMATClX, ’EVAPORATOR  TEMP  =’ ,F8. 2, 5X, ’CONDENSER  TEMP  =’,F8.2) 

41  FORMATClX, ’RADIATION  PARAMETER  =’,F8.2,  5X,’HEAT  TRANS  COEFF  =’, 

1  F18.2,  5X, ’PRESSURE  DROP  =’,F8.2) 

42  FORMATClX, ’TPE,  EKP,  EKL,  TPWE,  TVE, TEMPDIFF ,  BETAl,  UHPP,  PKPC’) 

43  FORMATC  1X,F8.2,  4(4X,F18.2)  ) 

44  FORMATClX, ’TV,  VDENS,  VISCOSV,  HFG,  TVC,  CKP,  CKL,  TPWC,  TPC’) 

45  FORMATClX,  4(4X,F8.2)  ) 

46  FORMATClX,  5F15.2) 
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ifHHUft*********************************** 

* 

*  PART  2  :  HEATPIPE  -  FIN  RADIATION 

* 

fH^i(0^:^^HHL**0***************>H************ 

* 

iL^fm^^^fUl^L^t****************************** 

*  CONSTANTS  RELATED  HEATPIPE-FIN 

RSTAR  =  RRO  /  HFLEN 
TSTAR  =  HTHICK  /  HFLEN 
PARAA  =  DSQRT(  0.4D0*BETA  ) 

RHOP  =  1.0  -  EPSIP 

RHOF  =  1.0  -  EPSIF 

PIPENO  =  NPIPE  -  1 

FINNO  =  (  NFIN  -  1  )  /  2 

NFINl  =  NFIN  -  1 

NFIN2  =  (  NFIN  +  1  )  /  2 

NDFIN  =  NFIN  /  8 

NDFIN2  =  NFIN  /  20 

NFIN3  =  NFIN2  -  NDFIN2 

NDPIPE  =  NPIPE  /  9 

ZETAD  =  l.ODO  /  FLOAT(FINNO) 

P  =  SQRTC  RSTAR**2  -  TSTAR**2  ) 

THETAMIN  *  ATAN(  TSTAR  /  P  ) 

THETAD  =  (  PI/2.0  -  THETAMIN  )  /  PIPENO 
THETAD2  =  0.5  *  RSTAR  ♦  THETAD 
THETAD3  =  THETAMIN  +  (  THETAD  /  4.75  ) 
ARCONST  =  FINNO  *  THETAD  ♦  RSTAR 
THETA  =  THETAMIN 

♦  . 

NCOUNT  =  1 

VSINE(l)  =  SIN(  THETAD3  ) 

VCOSIN(l)  =  COS(  THETAD3  ) 

DO  200  J=2,  NPIPE 

THETA  =  THETA  +  THETAD 
VSINE(J)  =  SIN(  TEHTA  ) 

VCOSIN(J)  =  COS(  THETA  ) 

200  CONTINUE 
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♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦***************** 

*  VIEW  FACTOR  FROM  FIN  SEGMENT  TO 

*  PIPE  SEGMENT  (  VFPF  ) 

iHc*if**rtiiti:tL:tHit**************************** 

* 

300  QJ3  =0.0 

* 

DO  329  J=l,  NPIPE 
IF  (  QJ3  )  326,  321,  321 

321  THICKTAN  =  (  TSTAR  *  VSINE(J)  )  /  VCOSIN(J) 

AJ3  =  RSTAR  ♦  VSINE(J)  -  TSTAR 
BJ31  =  THICKTAN  +  RSTAR  ♦  VCOSIN(J) 

BJ32  =  (  TSTAR  /  VCOSIN(J)  )*(  TSTAR  /  VCOSIN(J)  ) 

S  =  P  +  THICKTAN  +  (  ZETAD  /  5.0  ) 

AXJ3  =  S  *  VCOSIN(J)  -  TSTAR 
AXJ4  =  AJ3  *  AXJ3 

* 

IF  (  AXJ4  )  327,  328,  328 

327  VFPF(1,J)  =  0.0 
QJ3  =  -  1.0 

GO  TO  329 

328  B3  =  2.0  *  ((  S**2  -  2.0*S*BJ31  +  RSTAR**2  +  BJ32  )**  1.5  ) 
VFPF(1,J)  =  (  AXJ4  /  B3  )  *  ZETAD 

GO  TO  329 

326  VFPF(1,J)  =  0.0 

329  CONTINUE 

* 

QJ3  =  0.0 

* 

DO  311  J=l,  NPIPE 
QX3  =  0.0 

IF  (  QJ3  )  9999,  301,  309 

301  THICKTAN  =  (  TSTAR  *  VSINE(J)  )  /  VCOSIN(J) 

AJ3  =  RSTAR  *  VSINE(J)  -  TSTAR 

BJ31  =  THICKTAN  +  RSTAR*VCOSIN(J) 

BJ32  =  (  TSTAR  /  VCOSIN(J)  )  *  (  TSTAR  /  VCOSIN(J)  ) 

* 

DO  308  M=l,  NFINl 

IF  (  QX3  )  9999,  302,  307 

302  L  =  NFIN  +  1  -  M 
LX  =  L  -  1 

S  =  P  +  THICKTAN  +  LX*ZETAD 
AXJ3  =  S  ♦  VCOSIN(J)  -  TSTAR 
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AXJ4  =  AJ3  *  AXJ3 


IF  (  AXJ4  )  304,  303,  303 

303  B3  =  2.0  *  (  (  S**2  -  2.0*S*BJ31  +  RSTAR**2  +  BJ32  )**1.5  ) 
VFPF(L,J)  =  (  AXJ4  /  B3  )  *  ZETAD 

GO  TO  308 

t 

304  IF  (LX  -  2*FINN0  )  306,  305,  9999 

305  QJ3  =1.0 

306  QX3  =  1.0 

307  VFPF(L,J)  =  0.0 

308  CONTINUE 
GO  TO  311 

309  DO  320  L=l,  NFIN 
VFPF(L,J)  =  0.0 

320  CONTINUE 
311  CONTINUE 

SUHVFl  =0.0 

DO  313  N=2,  NFIN 
SUMVFl  =  VFPF(N,1)  +  SUMVFl 
313  CONTINUE 

VBETAl  =  ATAN(  (RSTAR*VSINE(1))  /  (2.0  +  P  -  RSTAR*VCOSIN(l) ) ) 
GAMMA  =  1.5708  -  VBETAl  -  THETAD3 

VFPF(1,1)  =  1.0  -  COS(GAMMA)  -  2.0*SUMVF1  +  VFPF(  NFIN,  1  ) 
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*  VIEW  FACTOR  FROM  HEAT  PIPE  TO 

*  HEAT  PIPE  (  VFPP  ) 

* 

* 

350  DO  357  J=1.NPIPE 
QK3B  *  0.0 

* 


DO 

357  K=l, 

NPIPE 

IF 

(  J  -  K  ) 

361, 

361,  362 

361  IF 

(  QK3B  ) 

9999, 

351,  356 

351  A351 

=  (  VSINE(J)  - 

VSINE(K)  )  *  RSTAR 

B351  =  ((1.0  +  P)*2.0)  -  (  RSTAR*(  VCOSIN(J)  +  VCOSIN(K)  )  ) 
E351  =  (  B351  ♦  VCOSIN(J)  -  A351*VSINE( J)  )*(  B351  *  VCOSIN(K) 
1  A351  *  VSINE(K)  ) 

* 

IF  (  E351  )  354,  353,  353 

353  C351  =  THETAD2  /  (  (  A351  **2  +  B351  **2  )♦*  1.5  ) 

VFPP(J,K)  *  E351  *  C351 

GO  TO  357 

354  QK3B  =1.0 

356  VFPP(J,K)  =  0.0 
GO  TO  357 

362  VFPP(J,K)  =  VFPP(K,J) 

357  CONTINUE 

DO  363  J=l,  NPIPE 
VFPP(J,1)  =  0.5  *  VFPP(J,1) 

363  CONTINUE 

* 

*  EDGE  TEMPERATURE  —  BLACK  SIMPLE  FIN 

* 

299  TEDGE  =  ((((  -2.2001614E-6  ♦  BETA  +  1.7217342E-4  )  *  BETA  - 

1  4.3974392E-3  )  ♦  BETA  +  4.6408689E-2  )  ♦  BETA  -  0.22322417  ) 

2  *  BETA  ♦  0.95864883 

* 

*  INITIAL  CHARACTERISTICS 

i^:^^itHHHHL^L0^Ll^^HL^iiHHLt********************* 

* 

TEMISMP  =0.0 
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100  DO  101  3=2,  NPIPE 

TEMISMP  =  VFPF(  NFIN2,  J  )  +  TEMISHP 
PRAD(J)  =  EPSIP 

101  CONTINUE 

TEMISMP  =  2. 0*ARC0NST*EPSIP* (TEMISMP  +  0.5*VFPF(NFIN2, 1) ) 
PRAD(l)  =  0.5  *  EPSIP 

IF  (  TEDGE**4  -  TEMISMP  )  103,  104,  104 

103  TEDGE  =  TEMISMP  *♦  0.25 

104  TEDGE5  =  TEDGE  **  5 
TEC  =  -  1.0 
1=1 

TF(I)  =  1.0 

106  IF  (  TEC  )  111,  112,  112 

111  TGRAD  =  (  TF(I)  **  5  )  -  TEDGE5 

IF  (  TGRAD  )  107,  107,  108 

107  GRADT(I)  =  0.0 
TEC  =  TF(I) 

GO  TO  109 

112  TF(I)  =  TEC 
GRADT(I)  =0.0 
GO  TO  109 

108  GRADT(I)  =  -  PARAA  *  SQRTC  TGRAD  ) 

TF(I+1)  =  TF(I)  +  GRADT(I)  *  ZETAD 

109  1=1+1 

113  IF  (I  -  NFIN2  )  106,  106,  399 

399  PARAB  =  (  RSTAR  *  EPSIP  )  /  (  1.0  +  P  ) 
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♦  IRRADIATION 
*************************^^****1^1********* 

* 

400  DO  402  L=l,  NFIN2 
FIRRAD(L)  =0.0 

M  =  NFIN  +  1  -  L 

* 

DO  401  J=l,  NPIPE 

FIRRAD(L)  =  PRAD(J)  ♦  (  VFPF(L,J)  +  VFPF(M,J)  )  +  FIRRAD(L) 

401  CONTINUE 

FIRRAD(L)  =  ARCONST  *  (  FIRRADCD  -  0.5*PRAD(1)  *  (  VFPF(L,1) 

1  +  VFPF(M.l))) 

IF  (  NCOUNT  -  10  )  410,  411,  411 

410  RF  =  NCOUNT 

FIRRADCD  =  RF  *  FIRRAD(L)  /  lO.ODO 

411  CONTINUE 

FIRRAD(M)  =  FIRRAD(L) 

402  CONTINUE 

* 

*  FIN  TEMPERATURE 

******^iJtit***************^t**********iti**** 

* 

500  N5  =  -  2 
511  TI  =  0.0 

* 

DO  501  L=l,  NFIN2 

TI  =  (  FIRRADCD  *  GRADTCD  )  +  TI 

501  CONTINUE 

506  SUMIRRAD  =0.0 
P3  =  0.0 

* 

DO  508  L=l,  NFIN2 
T4  =  TFCD  **  4 

GRADT2CD  =  BETA  *  C  T4  -  FIRRADCL)  ) 

522  SUMIRRAD  =  C  FIRRADCL)  ♦  GRADTCL)  )  +  SUMIRRAD 

GRAD  =  C  T4  *  TFCL))  -  TEDGE5  ♦  C  5.0  ♦  ZETAD+CTI  -  SUMIRRAD)) 

* 

IF  C  GRAD  )  504,  504,  505 

504  GRADTCD  =0.0 
P3  =  P3  +  1.0 
GO  TO  508 


12 


505  GRADT(L)  =  -  PARAA  *  SQRT(  GRAD  ) 

TF(L+1)  =  TF(L)  +  GRADT(L)  *  ZETAD  +  0.5*GRADT2a)*(  ZETADf+Z  ) 

508  CONTINUE 

CALL  COPPER(TPC,CKP) 

QFN  =-CKP*WC*HTHICK*(GRADT(l)/4,)*TPC/HFLEN 
N5  =  N5  +  1 

IF  (  N5  )  507,  511,  509 

507  TEDGE  =  0.5  *  (  TF(NFIN2)  +  TEDGE  ) 

TEDGE5  =  TEDGE  ♦♦  5 
GO  TO  511 

509  RATIO  =  GRADTC  NFIN2  )  /  GRADT(  1  ) 

PR  =  RATIO 

P5  =  TEDGE  -  TF(  NFIN2  ) 

RD  =  PR  -  RATIO 
NCOUNT  =  NCOUNT  +  1 
TEDGE  =  TF(  NFIN2  )  -  P5 
TEDGE5  =  TEDGE  ♦♦  5 
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7tf***^tllfjtilll}HL^HL*t************************ti* 

*  RADIOSITY 

^i*^L^Hfi^fi*^fi**l|L***^L^ir^i|i^^i^^:H^^L^^:t‘*************** 

* 

600  DO  601  L=l,  NFIN2 

FRAD(L)  =  EPSIF  *  (  TF(L)  ♦♦  4  )  +  RHOF  *  FIRRAD(L) 

M  =  NFIN  +  1  -  L 
FRAD(M)  =  FRAD(L) 

601  CONTINUE 
QTIRRAD  =0.0 

* 

DO  604  J=l,  NPIPE 
TIRRAD  =0.0 

* 

DO  602  L=l,  NFIN 

TIRRAD  =  FRAD(L)  *  VFPF(L,J)  +  TIRRAD 

602  CONTINUE 

TIRRAD  =  TIRRAD  -  0.5D0  *  FRAD(l)  *  (  VFPF(1,J)  +  VFPF(NFIN,J)  ) 

* 

T2IRRAD  =  O.ODO 
DO  603  K=l,  NPIPE 

T2IRRAD  =  PRAD(K)  *  VFPP(J,K)  +  T2IRRAD 

603  CONTINUE 

PRAD(J)  =  EPSIP  +  RHOP  *  (  TIRRAD  +  T2IRRAD  ) 

QTIRRAD  =  TIRRAD  +  T2IRRAD  +  QTIRRAD 

604  CONTINUE 

QC=QFN+SIGMA*TPC**4.*WC*RR0*(PI/2.-THETAMIN-(QTIRRAD/4.)) 

QTT1=SIGMA*TPC**4 . * (HFLEN+SQRT(RR0**2 . -HTHICK**2 . ) ) *WC 

ETAEFF1=QC/QTT1 

QDIV=ABS( ( (QREQ/4 . ) -QC) /QC) 

WRITE(6,*)  QDIV,QREQ,QC,ETAEFF1 
IF(QDIV.LE.lE-6)  GO  TO  899 
QREQ=4.*  QC 
GO  TO  799 

* 

**************************************** 

*  CRITEIA 


* 


IF 

(  NCOUNT  - 

15  ) 

400,  400,  704 

704 

IF 

(  NCOUNT  - 

100  ) 

705,  800,  800 

705 

IF 

(  ABS(  RATIO  )  - 

0.015  )  706,  800,  400 

706 

IF 

(  P3  -  1.0 

)  800,  800,  400 
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* 

*  EFFICIENCY 

* 

899  CALL  STEAMDEN(TV,VDENS) 

CALL  STEAMHFGCTV.HFG) 

CALL  COPPERCTV.EKP) 

CKP=EKP 

CALL  WATER (TV, EKL) 

EKW=EKP 

PKW= . 185*WF*EKW+WD*EKL 

PKEE= (WD*WF*EKL*EKW+W*EKL*PKW) / ( (W+WF) *PKW) 

CKW=CKP 

CKL=EKL 

PKEC= (W*CKL+WF*CKW) / (W+WF) 

DTT=CONSTA/(WE*EKP)+CONSTB/(WE*PKEE)+(TV*PDELTA)/(VDENS*HFG) 
♦  +CONSTB/ (WC*PKEC) +CONSTA/ (WC*CKP) 

UHP=qREQ/ (DTT*PI*RR0**2) 

WRITE(6,*)  DTT.UHP 
800  ETAFIN  =  -  GRADT(l)  /  BETAl 

ETAEFF  =  PARAB  ♦(1.5708D0  -  THETAMIN  -  (QTIRRAD*  THETAD)) 

1  +  (  ETAFIN  /  (  l.ODO  +  P  )) 

RN  =  RSTAR  *  (  BETAl  **  0.3333  ) 

COUNT  =  NCOUNT 
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*i^**t**************>lt******************** 

*  OUTPUTS 

* 


WRITE(6,952) 
WRITE(6,903) 
WRITE(6,911) 
WRITE(6,902) 
WRITE(6,912) 
WRITE (6, 902) 
WRITE(6,913) 
WRITE (6, 902) 
WRITE(6,914) 
WRITE(6.902) 
WRITE(6,915) 
WRITE (6, 902) 
WRITE(6.916) 
WRITE (6, 902) 
WRITE(6,903) 
WRITE(6,904) 


TEDGE.PS.GRADTCl) ,GRADT(NFIN2) .RATIO 

TSTAR , RSTAR , EPSIP , EPSIF , BETAl , ETAEFF , ETAFIN , COUNT 

(  TF(N),  N=l.  NFIN2,  NDFIN2  ) 

(  GRADT(N),  N=l,  NFIN2,  NDFIN2  ) 

(  GRADT2(N),  w=i,  NFIN2,NDFIN2  ) 

(  FIRRAD(N).  N=l,  NFIN2,  NDFIN2  ) 

(  FRAD(N),  N=l,  NFIN2,  NDFIN2  ) 

(  PRAD(J),  J=l,  NPIPE,  NDPIPE  ) 

RD,  TI,  TENS,  QTIRRAD,  FINNO,  PIPENO,  RN,  BETA 
P , ARCONST , THETAMIN .PARAB , ZETAD , THETAD , THETAD3 


901  FORMATC  IX,  10F12.8) 

902  FORMATC  IX,  11F10.6) 

903  FORMATC  IX,  4F14.8,2I5,2F14.8) 

904  FORMATC  IX,  7F14.8) 

911  FORMATC  IX,  ’TEMPERATURE’,  2X,  ’DISTRIBUTION’  ) 

912  FORMATC  IX,  ’ TEMPERATURE ’, 2X ,  ’GRADIENT’ ,2X,  ’DISTRIBUTION’) 

913  FORMATC  IX,  ’SECOND’,  2X,  ’DERIVATIVE’,  2X,’0F’,  2X, 

1  ’TEMPERATURE’  ) 

914  FORMATC  IX,  ’IRRADIATION’ ,2X, ’ON’ ,  2X,’FIN’  ) 

915  FORMATC  IX,  ’FIN’,  2X,  ’RADIOSTY’  ) 

916  FORMATC  IX,  ’TUBE’,  2X,  ’RADIOSTY’  ) 

952  FORMATC  IX,  ’EDGE  TEMP=’,  F9.6,  5X,  ’P=’,F9.6,  5X, 

1  ’ROOT  GRAD=’,  F9.6,  5X,  ’TIP  TEMP  GRAD=’,  F9.7, 

2  5X,  ’GRAD  RATIO=’,  PE10.3  ) 

981  FORMATC  IX,  6F14.8) 


CLOSE CUNIT=3) 
CL0SECUNIT=6) 
GO  TO  9999 


9999  STOP 
END 
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**4c4<*4>**  ***  If  >tc**  ****4[«:tr*:|c4<«>4<******«***4< 


*  SUBROUTINE  WATER  * 

*  * 

*  FIND  CONDUCTIVITY  OF  WATER  * 

♦  BY  USING  LEAST  SQUARES  METHOD  ♦ 


SUBROUTINE  WATER(T,WK) 

DOUBLE  PRECISION  WK 
REAL  WK 

AA  =  -0.3985D6 
BA  =  0.5282D4 

CA  =  -0.6389D1 

WK  =(AA  +  BA*T  +  CA*T*T)*lE-6 

RETURN 

END 


*  SUBROUTINE  COPPER  * 

*  * 

*  FIND  CONDUCTIVITY  OF  COPPER  * 

*  BY  USING  LEAST  SQUARES  METHOD  ♦ 


♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦if*********** 


SUBROUTINE  COPPER(T,CC) 

DOUBLE  PRECISION  CC 
REAL  CC 

AB  =  0.3799D3 

BB  =  -0.3780D-1 

CB  =  0.1023D5 

CC  =  AB  +  BB*T  +  CB  /  T 

RETURN 

END 
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*  SUBROUTINE  STEAMDEN  ♦ 

♦  ♦ 

*  FIND  DENSITY  OF  STEAM  * 

*  BY  USING  LEAST  SQUARES  METHOD  ♦ 


SUBROUTINE  STEAMDEN (T, DENSITY) 

DOUBLE  PRECISION  DENSITY, SV 
REAL  DENSITY.SV 

AC  =  0.3956D-6 

BC  =  0.9677D3 

SV  =  AC*T  +(  BC/T  ) 

DENSITY  =  1  /  SV 

RETURN 

END 


0^L**^Li^Hll|L:^L^L*^Hi^itt***************************** 

*  SUBROUTINE  STEAMHFG  * 

♦  * 


*  FIND  LATENT  HEAT  OF  WATER  VAPORIZATION* 

*  BY  USING  LEAST  SQUARES  METHOD  * 


SUBROUTINE  STEAMHFG (T.HEATLATENT) 

DOUBLE  PRECISION  HEATLATENT 
REAL  HEATLATENT 

AD  =  0.1022D4 

BD  =  0.9197D1 

CD  =  -0.1531D-1 

HEATLATENT  =  (AD  +  BD*T  +  CD*T*T)*1E3 

RETURN 

END 
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*l^**t:**^tifHHfrttlti:^H‘*jHf:^ii,i,t*********************** 


*  SUBROUTINE  STEAMUG  ♦ 

*  * 

*  FIND  VISCOSITY  OF  STEAM  * 

*  BY  USING  LEAST  SQUARES  METHOD  * 


************************************* 


SUBROUTINE  STEAMUG(T, VISCOSITY) 

DOUBLE  PRECISION  VISCOSITY 
REAL  VISCOSITY 


AE  =  -0.2S57D-7 
BE  =  0.4088D-4 

VISCOSITY  =  (T  /  (  AE*T  +  BE  ))*1E-12 

RETURN 

END 


******************************************* 
*  CASE  INPUT  DATA  * 

******************************************* 


RRO  = 

0.01 

RV  = 

0.01 

WE  = 

0.15 

WA  = 

0.05 

WC  = 

1.50 

NG  = 

50 

NHP  = 

10 

TS  = 

0.0 

TV  = 

300 

ETG  = 

0.95 

FLEN  = 

0.01 

THICK 

*  0.005 

EPSIP 

=  1.0 

EPSIF 

O 

II 

NFIN 

=  101 

NPIPE 

*  31 
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Appendix  B.  VIEW  FACTOR  EQUATIONS 
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Figure  B.l.  Geomety  Used  in  Fin-to-Heatpipe  View  Factor  Determination 


B.l  View  Factor  from  Fin  to  Heatpipe  [3] 

The  geometry  of  the  problem  for  evaluating  the  view  factor  from  fin  to  heatpipe 
is  shown  in  Fig. B.l  and  the  paramerters  are  expressed  in  the  non-dimensional  terms 
which  is  divided  by  the  half  of  fin  length,  L 

Geometric  identities  : 


=  +  = 

(B.l) 

i;  s; 

=  K'  -  +  C  +  i; 

(B.2) 

sj  +  s;  _  i; 

-  L-,  -  s; 

(B.3) 

S' 

sinxp  = 

(B.4) 

S*2  =  52*"  +  5;" 

(B.5) 

Combining  Eqs.(B.2),(B.3),  and  (B.5)  yields 

^  L\sin9  —  S4 
^  [52*2  +  53*2]  V2 


(B.6) 
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Combining  Eqs.(B.l),(B.3),  and  (B.6)  yields 


sin^p  =  — 

V 


where 

u  =  L\sin6  —  5^ 
V  =  [1*2  -  2L\{rlcose  + 
Combining  Eqs.(B.l),(B.2),  and  (B.3)  yields 


(B.7) 


L’l  =  [rf  -  +  C  +  i'tane 


(B.8) 


From  Eq.(B.7) 


where 


d{sinip)  = 


vdu  —  udv 


dv 


du  =  sind 
_  L\  —  [r*co$d  4- 

V 


The  view  factor  between  two  infintesimal  segments  is  given  by 
(see  Ref  8,  Eq.(36)) 

T?  _  1  d{sini') 

-  l—dT^ 


(B.9) 


(B.IO) 


Combining  Eqs.(B.l),  (B.9),  and  (B.IO)  and  noting  that  increasing  (  corresponds  to 
decreasing  tp  yields 


1  [r*sm^  —  t’][L\cos9  —  C] 

=  2 - ^373 - ‘‘i 


(B.ll) 


where 

Z  =  L-Y-  2L\[r;cose  4-  ftane]  + 

cosU 
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and 


+  C  +  ftane 


The  view  factor,  Fi_.f(l,l),  from  the  first  heatpipe  condenser  surface  segment 
(r’A0/2)  to  the  immediately  adjacent  fin  segment(AC/2)  cannot  be  evaluated  with 
Eq.(B.ll).  To  evaluate  this  view  factor,  the  first  haetpipe  condenser  segment  was 
assumed  to  be  located  at  a  distance  along  the  heatpipe  condenser  of  r*A0/4.75  away 
from  the  fin  and  the  first  fin  segment  was  assumed  to  be  located  at  a  distance  along 
the  fin  of  A(/5.0.  The  view  factor,  F,^_i(l,l),  is  then  defined  so  that  the  sum  of  all 
the  view  factors,  from  the  fin  segments  to  the  point  on  the  heatpipe  condenser  at  a 
distance  r*A^/4.75  from  the  fin,  is  equal  to  the  value  that  would  be  obtained  for  the 
view  factor  from  the  fin,  as  a  single  surface,  to  the  same  point. 

The  error  in  replacing  the  continuous  fin-heatpipe  condenser  surface  by  finite 
segments  was  evaluated  by  adding  all  the  view  factors  to  a  given  point  and  comparing 
the  result  to  the  analytical  expression  for  the  view  factor.  The  maximum  error  was 
approximately  1.5  percent  which  indicates  good  results  with  a  negligible  error  since 
there  are  100  fin  segments. 

B.2  View  Factor  from  Pipe  to  PipefS] 

The  geometry  of  the  problem  for  evaluating  the  view  factor  from  heatpipe  to 
heatpipe  is  shown  in  Fig. B.2  and  the  parameters  are  expressed  in  the  non-dimensional 
terms  which  was  divided  by  the  half  of  fin  length,  L  ,. 

Geometric  identities  : 

xP  =  di  +  6  (B.12) 

(2* 

sinS  =  —  (B.13) 

c* 
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cos6  = 


c* 


a*  =  rl[sin6i  —  sin02 


6*  =  2[1  +  {r:^  _  r^y/2]  _  r;[cose,  +  cose^] 


c*2  =  q*2  +  6*2 


From  Eq.(B,12) 


cosip  =  cosB\cos8  —  sin0\sin6 


From  Eq.(B.lO) 


C’  ^ 

^2-1  =  -cosip~d92 


From  Eq.(B,12),(B.13),  and  (B.14) 


dtp 

df2 


(B.14) 

(B.15) 

(B.16) 

(B.17) 

(B.18) 

(B.19) 

(B.20) 


B-5 


Performing  the  differentiation  of  Eq.(B.20)  becomes 


r.*2 


where 


and 


da' 

w,  = 


db' 
Ml  = 


Combining  Eq.(B.13),(B.14),  and  (B.18)  yields 


cost/)  = 


b'cosOi  —  a'sinOi 


Combining  Eq.(B.19),(B.21),  and  (B.22)  yields 


^2^1  = 


[6‘cos^i  —  a*sm0i][fe*cos02  ~  a’sm^a]  ^IdB 

[a*2  -f  6*2]3/2  ^ 


where 


a*  =  r*[sfn^i  —  sinB^ 


and 


6*  =  2[1  +  {r*^  -  -  r*[cos0i  +  cos^j] 


(B.21) 


(B.22) 


(B.23) 
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Appendix  C.  TABULATED  RESULTS 
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1 

Table  C.l.  Characteristics  of  Radiator  Parameter  vs 

Condenser  Temperature 

Tv 

Conductivity 

300.00 

401.8 

0.00038 

320.00 

400.4 

0.00046 

340.00 

396.3 

0.00056 

360.00 

393.8 

0,00067 

380.00 

391.6 

0.00079 

400.00 

389.5 

0.00093 

420.00 

387.5 

0.00108 

Table  C.2.  Characteristics  of  temperature  Difference  and  total  heat  transfer  rate 
vs.  Operating  Temperature 

Tv 

Qtotal 

Tpe  -  Tpc 

300.00 

166.71 

6.814 

320.00 

201.73 

7.187 

340.00 

241.76 

8.349 

360.00 

278.23 

9.693 

380.00 

338.60 

11.253 

400.00 

396.35 

13.067 

420.00 

460.98 

15.183 
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Table  C.3.  Characteristics  of  Heat  Pipe  vs.  Operating  Temperature 


Tv 

^HP,p 

VtU 

300.00 

54914.8 

0.9722 

320.00 

57177.4 

0.9088 

340.00 

58992.6 

0.8546 

360.00 

60364.8 

0.8076 

380.00 

61297.2 

0.7672 

400.00 

61792.2 

0.7315 

420.00 

61851.5 

0.7000 

Table  C.4.  Characteristics  of  Black  Surface  Heatpipe  Radiators 

with  Fins  of  Neg- 

ligible 

Thickness 

• 

II 

o 

o 

o 

Cp  =  €/  =  1.0 

II 

o 

Beta 

midpt  Temp 

Efficiency 

Ef f iciency 

(fin) 

(eff) 

10.0 

0.636 

0.130 

0.723 

5.0 

0.676 

0.190 

0.747 

2 . 0 

0.751 

0.302 

0.793 

1.0 

0.815 

0.408 

0.836 

0.5 

0.874 

0.522 

0.882 

II 

o 

10.0 

0.580 

0.145 

0.600 

2 . 0 

0.726 

0.332 

0.702 

0.5 

0.864 

0.567 

0.832 

t,*  =  0.4 

10.0 

0.545 

0.155 

0.505 

2 . 0 

0.712 

0.352 

0.633 

0.5 

0.858 

0.594 

0.793 
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Appendix  D.  SIMPLE  FIN  SOLUTION 


D-l 


The  characteristics  of  a  simple  fin  with  no  incident  irradiation  were 
computed  using 


(D.l) 


The  fin  efficiency  is  given  by 

- 

0  dx* 


(D.2) 


Equation(D-l)  was  solved  by  iteration.  The  iterations  were  terminated  when  the  two 
computed  edge  temperatures(  Te)  were  within  5  x  10“^.  The  dimensionless  length 
of  the  fin  segment  was  equal  to  0.005.  The  results  are  given  in  Table  D.l.  The 
relation  between  the  edge  temperature(  Te  )  and  the  radiation  parameter(/?  )  can 
be  expressed  by  the  following  poiyiiomiai; 


Te  =  0.959  -  0.223/3  +  4.641  x  10_2(/3)^  -  14.347  x  10"^(/3)^ 


+  1.722  X  -  2.200  x  10-®(/3)* 


(D.3) 
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Table  D.l. 


Beta 


40,0 
20.0 
10.0 
4.0 
2.0 
1.0 
0.40 
0.20 
0.10 
0 . 04 


Charateristics  of  Simple  Fin  Radiators 


midpt  Temp 


0.335 

0.407 

0,487 

0.603 

0.693 

0.778 

0.871 

0.923 

0.956 

0,981 


Efficiency 


0.010 

0.141 

0.197 

0.303 

0.410 

0.535 

0.706 

0.815 

0.895 

0.958 
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Appendix  E.  PROPERTIES 


This  appendix  shows  the  experimented  data  selected  from  reference  18  and  the 
estimated  data  calculated  by  the  least  square  method.  The  error  of  all  properties 
between  the  experimened  data  and  the  estimated  data  lies  within  maximum  1.7  9c 
in  the  operating  temperature  rang  (  200  -  500  K  ). 
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Figure  E.l.  Thermal  Conductivity  of  Pure  Copper 
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Figure  E.2.  Density  of  Saturated  Water  (  Steam  ) 
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Figure  E.3.  Latent  Heat  of  Vaporization  (  Water  ) 
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Figure  E.4.  Viscosity  of  Water  (  Steam  ) 


E-5 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

lO 

o 

in 

o 

in 

o 

m 

o 

m 

o 

to 

to 

m 

m 

'g- 

'g 

cn 

cn 

(N 

CM 

Aimonpuoo  /eu/jet/1 


Figure  E.5.  Thermal  Conductivity  of  Water  (  Liquid  ) 
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